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Abstract 

We use molecular dynamics computer simulations to study the relaxation 
dynamics of Na20-2(Si02) in its molten, highly viscous state. We find that 
at low temperatures the incoherent intermediate scattering function for Na 
relaxes about 100 times faster than the one of the Si and O atoms. In contrast 
to this all coherent functions relax on the same time scale if the wave-vector 
is around lA _1 . This anomalous relaxation dynamics is traced back to the 
channel-like structure for the Na atoms that have been found for this system. 
We find that the relaxation dynamics for Si and O as well as the time depen- 
dence of the coherent functions for Na can be rationalized well by means of 
mode-coupling theory. In particular we show that the diffusion constants as 
well as the a-relaxation times follow the power-law predicted by the theory 
and that in the /3-relaxation regime the correlators obey the factorization 
property with a master curve that is described well by a von Schweidler-law. 
The value of the von Schweidler exponent b is compatible with the one found 
for the mentioned power-law of the relaxation times/diffusion constants. Fi- 
nally we study the wave-vector dependence of f s (q) and /(<?), the coherent 
and incoherent non-ergodicity parameters. For the Si and O atoms these 
functions look qualitatively similar to the ones found in simple liquids or pure 
silica, in that the coherent function oscillates (in phase with the static struc- 
ture factor) around the incoherent one and in that the latter is approximated 
well by a Gaussian function. In contrast to this, f(q) for Na-Na is always 
smaller than f s (q) for Na and the latter can be approximated by a Gaussian 
only for relatively large q. 
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I. INTRODUCTION 



In the last two decades our understanding of the structural and dynamical properties of 
glass forming liquids has increased impressively . This progress is due to significant ad- 
vances in various experimental techniques (light and neutron-scattering, dielectric measure- 
ments, etc.), the development of new theoretical approaches and concepts (mode-coupling 
theory (MCT), landscapes, etc.) and last not least to the extraordinary advances that 
computer simulations have made PHT3f. The results of all these efforts is a widely accepted 
picture on the relaxation dynamics of glass-forming liquids: At high temperatures this dy- 
namics is the one of a normal liquid and hence the typical time correlation functions decay 
exponentially in time and the T-dependence of the typical relaxation times r(T) shows an 
Arrhenius-law. (Note that although there is a well developed theoretical machinery to de- 



scribe the static properties of such liquids [14], an accurate understanding of the dynamics 



is still lacking.) With decreasing temperature the dynamics changes in that time correlation 
functions are no longer exponential and the T-dependence of the relaxation times is often 
super-Arrhenius. In the past it has been shown that the slow dynamics observed in this 
temperature range can be rationalized very well by means of the so-called mode-coupling 
theory of the glass transition, and that often this theory is not only able to describe this 
dynamics qualitatively but even quantitatively [|,|T^-f2T]] . If the temperature is decreased 
even further the relaxation times increase very quickly and typically show an Arrhenius 
dependence on T with an activation energy that is larger than the one found at higher tem- 
peratures. The details of the dynamics in this temperature range are not understood very 
well and also there exist no satisfactory theoretical description for it. 

The presence of the three mentioned regimes in r(T) (Arrhenius, super-Arrhenius, Ar- 
rhenius) is the typical behavior found in the so-called "fragile" glass-forming liquids J22 



which include polymers, most molecular glass-formers, and simple liquids. In contrast to 
this, the so-called "strong" glass-forming liquids (typical examples include silica and many 
other oxide glasses) show in the whole accessible temperature range an Arrhenius-law, i.e. 
no super-Arrhenius temperature dependence is found. Since it is one of the main achieve- 
ments of MCT to rationalize this super-Arrhenius T-dependence, it was believed for quite 
some time that the theory is not very useful for these type of glass-formers. In recent com- 
puter simulations it was shown, however, that most models for strong glass-forming liquids 
exhibit at sufficiently high temperatures deviations from the Arrhenius law found at inter- 
mediate and low temperatures [p3-28[ and that in the temperature range in which these 



deviations are seen the relaxation dynamics can be described very well by MCT [25,29,30 



in agreement with the results of some experiments |3l|j32][ . Therefore one can conclude that 
this theory is not only able to describe fragile glass-formers, but also all the intermediate 
and strong ones. 

The results discussed so far concern systems in which the dynamics of all atomic species 
occurs on roughly the same time scale. There are, however, large classes of materials in 
which this is not the case. E.g. in superionic glass-formers, such as Na20-xSi02 or mixed 
alkali glasses such as 0.5Na2O-0.5K 2 O-3SiO2, the alkali atoms move on a time scale which 
at low temperatures is many orders of magnitude faster than the one of the atoms consti- 
tuting the matrix (here Si and O) |33|-|3"8[| . To what extent MCT is able to rationalize the 
very heterogeneous dynamics of these type of glass-forming systems is presently unknown, 
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although recently evidence has been found that certain aspects of the relaxation dynamics 
can indeed be understood by means of the theory [39]. The goal of the present paper is 
therefore to do a detailed investigation of the relaxation dynamics of a prototype of such 
an ion-conducting glass-former, Na 2 0-2Si0 2 , and to see to what extent the predictions of 
MCT regarding the dynamics of glass-forming liquids hold. Note that although so far there 
has been no test of MCT for these systems, computer simulations have already been used for 
a long time to study such ion-conducting materials. Seminal work on this goes back more 
than twenty years to Soules et al. |40| who studied the structure of sodium-disilicate and to 
Angell et al. who looked at certain aspects of the diffusion dynamics f[|T],[EJ. These early 
investigations were more recently followed up by similar studies on related systems and also 
many experimental investigations |43|-|63|. 

A further interesting feature of these systems is that the distribution of the alkali ions 
in the Si02 matrix is not completely uniform but instead forms an interconnected network 

m 



of small pockets and filaments 



Although for a long time there was no 



direct experimental evidence for the presence of these structures, recent neutron-scattering 
experiments of Meyer et al. on Na20-2Si02 showed a structural feature at a wave-vector 
q p& 0.9 A -1 [Q, in excellent agreement with the results from computer simulations that 
could attribute this peak to the presence of the above-mentioned network |]3Tp3^ , |B"6| . How 
the presence of these channels affects the relaxation dynamics is, however, so far not known 
and therefore the present paper is also a contribution to gain insight into this matter. 

The remaining of the paper is organized as follows: In the next section we will describe 
the model used for the simulations as well as give their details. The following section is then 
devoted to the presentation of the results and in the final section we summarize and discuss 
them. 



II. MODEL AND DETAILS OF THE SIMULATION 



The force field used in the present work is a modification of the one proposed by Kramer et 
to describe zeolites |)7|], i.e. crystalline materials that contain also Si, O and Na. This 



al. 

potential had the functional form 



<t>*p{r) 



q a q/3e' 
r 



a 



+ A aP exp (-B a pr) - a, f3 e [Si, Na, O] 



(1) 



where r is the distance between two atoms of type a and (3. The parameters A a p, B a p, and 
C a p can be found in Ref. |67| . Note that q a is an effective charge with values gsi = 2.4 and 
go = —1.2. In Ref. [67| the effective charge of the sodium atoms was chosen to be q^a, = 1.0, 
which has the disturbing effect that Na 2 is not neutral. Hence we changed the potential 
of Ref. |67| somewhat, by assigning to the Na atom an effective charge q^ & = 0.6 p6|. This 



change of the charge affects of course the local structure and therefore we have added to 
<p a p(r) a term which compensates this change at short distances. Hence the potential we 
used in our simulation is given by 



ap\ 



<Mr) + [1 



(1 - <W)(1 - Sp Na )} 6(r c - r) 



(2) 
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with gsi = 2.4, go = —1-2, and g~Na = 0.6 In [C{r c — r) 2 + 1]. (Here 6 is the usual Heaviside 
function.) The parameters C = 0.0926 A~ 2 and r c = 4.9 A were chosen such that at 
ambient pressure the resulting structure is in good agreement with the experimental results 



from neutron scattering | 66|| . Thus the potential $ a ^(r) used in the present work is at short 
distances very similar to the one proposed by Kramer et al. and is at large distances (where 
the effective charges are screened) modified such that all systems of the form Na20-xSi02 
are neutral. 

In previous studies it has been shown that this potential is able to reproduce many struc- 
tural properties of Na20-:rSi02, with x = 2,3,4, and also certain aspects of the relaxation 
dynamics of this model are in good qualitative agreement with experiments p9] , |61| , |65| , |66| , |68| . 
Although it cannot be expected that such a simple model is able to reproduce faithfully 
all the features of the relaxation dynamics of the real material it can be expected that the 
salient properties are at least correct from a qualitative point of view. 

In the present simulation we integrated the equations of motions using the velocity 
form of the Verlet algorithm with a time step of 1.6 fs. At each temperature we first 
equilibrated the system by coupling it to a stochastic heat bath. This equilibration time 
exceeded the structural relaxation time of the silicon atoms (measured by means of the 
incoherent intermediate scattering function for a wave-vector 1.7 A -1 , which corresponds 
to the length scale of the distance between two tetrahedra, see below). Subsequently we 
started a production run in the NVE ensemble using a density of 2.37 g/cm 3 which is close 



to the experimental value |69|]. The number of particles was 8064 (As; = 1792, No = 4480, 
-^Na — 1792) in a cubic simulation box of size L = 48.653 A. Such a large system size is 
needed if one wants to avoid finite size effects in the relaxation dynamics JT0,7I|. In addition, 



previous simulations have shown that sodium-silicate glasses have structural features that 
extend over 6-8 A and therefore also this calls for rather large system sizes. The 

temperatures investigated were 4000 K, 3400 K, 3000 K, 2750 K, 2500 K, 2300 K, and 
2100 K. At the lowest temperature the length of the production run was 1.5 million time 
steps, which corresponds to 2.5 ns. In order to improve the statistics of the results we did 
at each temperature two completely independent runs. 



III. RESULTS 

The main quantities of interest in the present paper are -F(q, t) and -F s (q, t), the coherent 
and incoherent intermediate scattering functions for wave-vector q. (Note that for isotropic 
systems, like the one studied here, these space-time correlation functions depend only on q, 
the modulus of q, i.e. there is no directional dependence. Therefore we have made use of 
this fact also in the analysis of our data and have averaged over all wave-vectors with the 



same modulus.) These two observables are not only of great theoretical interest [1J] but can 



also be measured directly in neutron-scattering experiments |[72|| . They are given by 

F aP (q,t) = = £ £<exp [iq • (r,(t) - r fc (0))]> (3) 
iV j=ik=i 

i N a 

KM = - ]T(exp [zq ■ ( rj (t) - r,(0))]> . (4) 

iV 3=1 
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Here rj(t) is the position of particle j at time t, N a is the number of atoms of type a, and 
N is the total number of atoms. 

In Ref. J6£| we have shown that in NS2 the partial structure factors SgiSiiq), Sgiof^), 
and Soo(q) have two pre-peaks: one at q\ = 0.94 A -1 and a second one at g 2 = 1-7 A -1 . 
(The "main" -peaks, corresponding to the length scale of a nearest neighbor pair of Si-0 
and Na-0 are at m 2.8 A -1 and w 2.1 A -1 , respectively.) The peak at qi is related to the 
distance between neighboring tetrahedra and thus corresponds to the so-called "first sharp 
diffraction peak" in pure silica. The peak at q\ has recently been shown to be related to the 
typical distance between the channels mentioned in the introduction P9"|. Unfortunately, 



due to cancellation effects in the partial structure factors (weighted with the appropriate 
experimental neutron scattering cross sections) it is hard to see this peak in a neutron 
scattering experiment done at room temperature [E3|JjH,[73] . However, the recent neutron- 
scattering studies of Meyer et al. have shown that this system shows at high temperatures, 
i.e. above T « 1200 K which is well above the glass transition temperature p3, a feature 



in the elastic signal at a wave-vector around 0.9 A" 1 , i.e. very close to q\ [64]. Thus we 
conclude that this structural feature is not only seen in the present model for NS2, but 
can be found in the real material as well. In the following we will demonstrate that these 
structural features have also a strong influence on the relaxation dynamics of the system. 

In Fig. ^ we show the time dependence of the incoherent intermediate scattering function 
for all temperatures investigated. We see that at high temperatures the relaxation of the 
structure is relatively fast in that the correlation function decays, after the microscopic time 
scale which lasts around 0.2 ps, basically exponentially. From this graph we also recognize 
that at these temperatures this decay is about a factor of four faster for the sodium atoms 
than the one for the silicon atoms, and about a factor of three faster than the one for the 
oxygen atoms. This is in agreement with the values for the diffusion constants which at 
this temperature show the same dependence on the species f66 |. At low temperatures the 
correlators have a time dependence that differs qualitatively from the one at high T in that 
they show a plateau at intermediate times. This plateau is related to the fact that on 
this time scale the particles are trapped by their surrounding neighbors (a behavior which 
is often called "cage-effect") and hence the time correlation function changes only slowly. 
Only for relatively large times the particles are able to leave this (temporary) cage and 
hence the correlators finally decay to zero. It is customary to call the dynamics in the 
time window in which the correlators are close to the plateau the "/^-relaxation" whereas 
the second relaxation step is called the "a-relaxation" . From the figure we see that at low 
temperature, r(q) (the time scale for the a-relaxation) of the sodium atoms is about a factor 
of 30 smaller than the one for the oxygen atoms and about a factor of 100 smaller than the 
one for silicon. Hence we conclude that the temperature dependence of r depends on the 
species, also this in agreement with the one found for the diffusion constant ^6f. Below we 
will discuss this dependence in more detail. 

Also included in the figure is a fit to the curves at the lowest temperatures with a 
Kohlrausch-Williams-Watts-law (KWW), i.e. F"(q,t) = v4exp(— {t/r) 13 ) (dashed lines). 
We see that this functional form is able to describe the data very well in the a-relaxation 
regime, as it is usually the case for the relaxation dynamics of glass-forming liquids. The 
value of the exponent (3 is 0.8 for Si and O and 0.47 for Na. Thus we see that the relaxation 
dynamics is not very stretched for the atoms making up the matrix whereas it is very 
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stretched for the network modifier. Below we will come back to a possible interpretation 
for this low value of (3. (Note that a high value of (3 for Si and O is in agreement with 
the observation that also in pure silica the correlators for these length scales are not very 
stretched [|S6|| . The presence of sodium does, however, indeed lower the value of (3 somewhat.) 



Having discussed the species and temperature dependence of F"(q, t) we now focus on the 
coherent function F al3 (q,t), defined in Eq. (^|), since that one gives the information how the 
overall structure of the system relaxes. The time and temperature dependence of F SlSl (q,t) 
and ^ NaNa (g ; t) is shown in Fig. |2] (the one for the 0-0 correlation is similar to the one for 
Si-Si). A comparison of the curves in this figure with the incoherent functions, Fig. [TJ, shows 
that for the case of silicon the coherent and incoherent functions are very similar in that 
the height of the plateau at intermediate times as well as the typical relaxation times for 
the a-relaxation are quite comparable. This is the usual behavior found for glass forming 



systems such as pure silica, water, or simple liquids [p0|,[74i|r5| . In contrast to this the time 
dependence of _F s Na (g, t) for sodium differs strongly for the coherent function F NaNa (g, t) in 
that the relaxation time of the former is about a factor of ten smaller than the one of the 
latter. This shows that the relaxation dynamics of the sodium atoms has unusual features. In 



fact it has been shown in Refs. [^,0] that the motion of a tagged Na atom is relatively fast 
in that the particles hop between certain preferential sites (generated by the Si-0 matrix). 
However, the relative spatial arrangement of the Na particles (i.e. their global structure) 
does not change under this type of motion and hence the coherent function decays much 
slower than the incoherent one. 

One of the important predictions of MCT is that close to the critical temperature of the 
theory the shape of the time correlation functions does not depend on temperature. This 
means that a correlator 0(t, T) can be written as 

0(t,T) = 0(t/r(T)) , (5) 

where r(T) is the a-relaxation time at temperature T. To what extent this prediction, 
which is often called the time-temperature superposition principle (TTSP), holds can easily 
be tested by plotting the correlators versus t/r. For this we have defined the a-relaxation 
time r as the time it takes the correlator to decay to 0.1. Although this definition (and 
the value 0.1) is somewhat arbitrary, it is easy to see that if the TTSP holds the details 
of the definition do not matter. In Fig. |3| we show the correlators as a function of t/r for 
all temperatures investigated. We see that for silicon and oxygen the TTSP holds basically 
for the whole T-range. In the case of sodium we have already shown that its dynamics is 
very different from that of Si and O in that, e.g., for F s Na (g, t) one obtains a very different 
stretching exponent (3 from that for Ff 1 and F® at a given q. However also for F s Na (g, t) the 
TTSP work for low temperatures, T < 2500 K. That for this species the TTSP does not 
work for higher temperatures is reasonable since for T > 2750 K the typical a-relaxation 
times are still of the order of a few picoseconds only (see Fig. |l]b) and thus the a-relaxation 
regime is not well separated from the microscopic dynamics. Note that if one would be able 
to equilibrate the system at even lower temperatures it can be expected that the TTSP 
starts to break down, since at low T the dynamics is likely to be dominated by a simple 
diffusive motion in which the single step is the breaking of a Si-0 bond, a behavior which 
has been found in pure silica [J25|J . Hence one can expect that this affects also the sodium 
dynamics and thus to the breaking down of the TTSP. 
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The correlators discussed so far were for q = 1.7 A" 1 , the length scale corresponding 
to the typical distance between two neighboring tetrahedra. We now investigate how the 
a-relaxation time depends on the wave-vector as well as on the particle species. In Fig. |] 
we show the g-dependence of r as determined from the incoherent correlators as well as 
some of the coherent ones. The temperature is T = 2100 K, thus the lowest temperature at 
which we were able to fully equilibrate the system. Since for the case of a diffusive motion 
the relaxation time r(g) of the incoherent function is equal to (Dq 2 )^ 1 (D: diffusion con- 
stant) we plot directly r(g)g 2 which can be understood as an inverse g-dependent diffusion 
constant. We first start with the particles forming the matrix, i.e. Si and O. We see that for 
these elements r(q)q 2 from the incoherent functions increases weakly with increasing q, goes 
through a maximum and then decreases again at large q. The location of this maximum is 
around 2.8 A -1 , i.e. at the wave-vector at which the static structure factor has the peak 
that corresponds to the distance between neighboring Si-O. The same type of g-dependence 
has also been found in the case of pure SiC>2 [ 30 1 and it reflects the fact that on this length 
scale the system is particularly rigid. The relaxation times of the coherent functions have a 
more complex g-dependence in that they oscillate around the r for the incoherent function. 
This oscillation is in phase with the corresponding static structure factor, an effect that is 
known as de Gennes-narrowing |L4] and which is found in simple liquids J75| but also in 



a system like pure silica ||30|| , or water |74|. From the figure we recognize that the largest 



relaxation time is at a wave-vector q « 0.95 A 1 , which corresponds to the length scale of 



the network of channels discussed in Refs. P9| , pl| . The different chemical ordering of silicon 



and sodium leads to the presence of an additional intermediate length scale which is re- 
flected in the dynamics of sodium in that the trajectories of the sodium atoms are restricted 
to a network of channels in a Si-0 matrix. This network of channels is reflected in the 
collective correlations by the slowest relaxation process of our system, the rearrangement 
of the channel structure. Finally we remark that for wave-vectors below 0.95 A -1 r(q)q 2 
exhibits a very steep increase: We see that in this range of q values the curves are very 
well compatible with a straight line which corresponds to a growth of the relaxation times 
like r(g) oc q~ 2 exp(Aq), where A is a positive constant. Although we are not aware of any 
theoretical reason for such a dependence it seems to describe our data remarkably well over 
two decades in r. 

For the sodium atoms the g-dependence of r is more complicated than the one for Si 
and O. For the coherent correlators this function is qualitatively similar to the one found 
for Si-Si and O-O. In particular we find again a pronounced peak at q ~ 0.95 A -1 , i.e. the 
length scale of the channels. Note that close to this peak also the absolute value of r is close 
to the ones for Si and O which shows that on this length scale the spatial arrangement of the 
Na atoms can only relax if the Si-0 matrix relaxes. The relaxation time for the incoherent 
Na function behaves very differently. First of all we see that it is significantly smaller than 
the one of the coherent function, in agreement with our conclusions from Figs. [l| and Also 
this function increases for small and intermediate wave-vectors, shows a maximum at the 
location of the corresponding peak in the partial structure factor (included in the figure as 
well) and then decreases for even larger q. It is not possible to determine the relaxation 
times of F NaNa (g, t) for g > 2.5 A -1 with our definition F af3 (q,t = r) = 0.1 because for 
large wave- vectors F NaNa (g, t) has already decayed to values around or lower 0.1 before the 
a-relaxation starts to develop (see Fig. || below). In contrast to that we can determine the 
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relaxation times for F^ a (q,t) also for relatively large wave-vectors. Since for sufficiently 
large q one has F*(q,t) ~ F aa (q,t), we know that in this limit the relaxation times for 
F aa (q,t) have to approach the ones for F°(q,t). This implies that the curve r{q) for the 
Na-Na correlation must start to decrease quickly in order to approach the one for the Na 
correlation. Hence we expect that for large wave-vectors the relaxation times for the Na-Na 
correlation become significantly smaller than the one of the matrix. An inspection of the 
correlation functions shows however, that this is not yet the case for wave-vectors smaller 
than 3.5 A -1 . 

In an earlier study of pure silica, the paradigm of a "strong" glass-former, we found 
the surprising result that the temperature dependence of the diffusion constant as well as 
the viscosity show at high temperatures a significant deviation from the expected Arrhenius 
law [p5|1 . In that paper it was argued that this deviation can be rationalized by mode- 
coupling theory which predicts the existence of a "critical temperature" T c close to which 
the transport coefficients show a non-Arrhenius behavior ||. Since significant deviations 
from an Arrhenius law have been observed also for the sodium-silicate system investigated 
here it is reasonable to see to what extent MCT is able to rationalize the relaxation 
dynamics. The theory predicts that close to T c the temperature dependence of the diffusion 
constants or a-relaxation times is given by a power-law: 

D(T) oc (T — T c ) 7 and r(T) oc (T — T c )~ 7 . (6) 

Here 7 is a system-universal constant, i.e. it does not depend on what species or wave-vector 
one considers. In addition MCT predicts that the value of the exponent 7 has a one-to-one 
correspondence with the exponent b of the so-called von Schweidler-law that is discussed 
below. This connection is given by 

o, 1 4. 1 with [ r ( x - a )] 2 F( 1 + 2& )] 2 m 

1 = Ya + 2b Wlth r(l-2a) = r(l + 26) ' (?) 

i.e. the second equation can be used to determine the value of a from b and then the first 
equation can be used to calculate 7. (Here T(x) is the usual T— function.) In Ref. |J9] we 
have shown that the value of b is around 0.47 (see also Fig. |6] below). Hence Eqs. (|7|) give 
a value of 7 = 2.87, which can be considered as the theoretical estimate of MCT for the 
exponent. If the prediction of the theory on the power-law and the value of the exponent is 
correct, a plot of t- 1 / 7 (or D 1 /" 1 ) vs. T should give a straight line. This type of plot is shown 
in Fig. |5| for the relaxation times r(q) for the wave-vectors q = 0.94 A -1 and q = 1.7 A -1 
as well as the diffusion constants (which were determined from the long time limit of the 
mean squared displacement of a tagged particle |66fl). Note that we have not included the 
data for Na, since, as shown in Ref. |)6] , the diffusion constant for Na follows an Arrhenius 



law in the whole temperature range, and hence the power-laws given by Eq. (y) certainly 
do not hold for this species. We see that all the curves do indeed show a straight line in 
a temperature interval that is quite substantial. Linear fits in this region are included in 
the graph as well (solid lines). Furthermore we see that the extrapolation of these straight 
lines to lower temperatures intersects the T-axis at a point, the critical temperature T c , 
which depends only weakly on the time considered, and which is around 2000 K. Hence 
we conclude from this figure that the relaxation times and diffusion constants for Si and 
O do indeed show the predicted power-law dependence with a common exponent 7 and a 
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common critical temperature T c . Finally we note that in the temperature interval in which 
t and D show this power-law these quantities change by about two orders of magnitude and 
hence the existence of the power-laws is not just in a trivial matter. Note that although 
MCT predicts that at T c the relaxation times should diverge, in reality this is not found. 
The reason for this is that once the relaxation times have increased beyond a certain value, 
for atomic systems usually on the order of 10 ns, the system starts to relax via processes 
that currently can be taken into account by the theory only in a schematic way. Despite 
the presence of these processes, usually called "hopping processes" the theory is still able 
to make prediction on the relaxation dynamics on the time scale of the /^-relaxation. For 
more details we refer to Refs. P,|76| . The presence of these hopping processes for the present 



system can be inferred from Fig. [| in that the data points for T < 2600 K are above the 
theoretical straight line, i.e. the relaxation is faster than predicted from the power-law. 

Having discussed the relaxation dynamics of the system in the a-regime, we now turn 
our attention to the /3-regime, i.e. the time window in which the correlators are close to the 
plateau (see Fig. |l|). One of the main predictions of MCT is that in the /3-regime and for 
temperatures close to T c the time dependence of the correlation function is system universal 
in that any time correlation function <fi(t) can be written as 

<f ) (t) = <f )c + h^G(t) . (8) 

Here <p c is the height of the plateau, often also called non-ergodicity parameter, and the 
whole time dependence is in the system universal function G(t). Due to its structure, 
Eq. (H|) is often called "factorization property", since 4>{t) — C factors into a time dependent 
function and a ^-dependent function. (Note that this factorization property holds also for 
the case that hopping processes are present P,|76|.) The time dependence of G{t) is given 
by the solution of a non-linear equation which can be solved numerically [H. However, it 



can be shown that very close to T c this solution is given by the following form [77] : 

G(t) = -Bt h + B't 2b . (9) 

This approximation is good for times at which the correlators have started to fall below 
the plateau, but are still close to it. The first power-law of Eq. (|9]) is often called "von 
Schweidler-law" and the exponent b the "von Schweidler-exponent" . Note that, since G(t) 
is predicted to be independent of the correlator, also the value of b should be the same for all 
4>(t). However, since Eq. (|9]) also includes the first correction term to this asymptotic result, 
the coefficient B' will depend on the correlator. In Fig. |B| we show the time dependence 
of the incoherent intermediate scattering function for the oxygen atoms for various wave- 
vectors (at T = 2100 K). Also included are fits with the functional form given by Eqs. (§) 
and ([]) using the height of the plateau as a fit parameter. The dotted lines correspond to 
a fit in which G(t) is only given by the von Schweidler-law, i.e. the first term on the RHS 
of Eq. (H), whereas the dashed line is the case that also the second term in Eq. (|9|) is taken 
into account. In these fits B and B' were fit-parameters that were allowed to depend on 
q, whereas the exponent b was a global fit-parameter. From the figure we conclude that 
the von Schweidler-law is indeed able to describe well the dynamics close to the plateau. 
Furthermore we see that the inclusion of the correction term increases the time window for 
which this law holds by about a factor of ten in qualitative agreement with the results of 
such an analysis for other glass-forming liquids p|, |30| , |7^ , [75| . We also mention that a similar 
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good fit is obtained for the case of F s (q,t) for Si. In Ref. |39| we showed the same type 
of fits for the coherent functions of Si and O and found that also for these correlators the 
/3-regime is described very well by the functional form given in Eq. (^). Hence we conclude 
that in the /3-regime the relaxation dynamics is indeed independent of the correlator. 

We emphasize that this universality holds only for the /^-relaxation regime and not for 
the a-relaxation. For the latter one finds, that the stretching exponent /3, and therefore the 
shape of the correlator, depends on the species or the wave-vector. In particular we have 



discussed in Ref. |39| the wave-vector dependence of (3 for the Na atoms and have shown 
that for small and intermediate g, q < 1.5 A -1 , (3 changes significantly thus showing that 
the a-relaxation is indeed not universal from this point of view. 

In this context it is, however, very interesting that (3(q) for the Na atoms becomes 
independent of q for q > 1.6 A -1 [j59| , although in this g-range the structure factor SWaNaC?) 
still shows pronounced features, i.e. is not a constant. Such a behavior was predicted some 
time ago by Fuchs who used MCT to show that lim g ^ 00/ 9(g) = b, i.e. for large wave-vectors 
(3 should converge to the von Schweidler exponent b |78[]. We have found that in our case 
this is indeed the case, i.e. that for large q the stretching exponent is indeed compatible 



with b = 0.47 fl39|. Hence this nice agreement between our results and the prediction of 
MCT shows that the theory is indeed able to describe also this features of the relaxation 
dynamics of the present system. We also mention that for the case of a hard sphere system 
the MCT prediction for the independence of (3 of q holds only for wave-vectors that are 
much larger than the location of the first peak in the structure factor. The fact that for the 
present system this asymptotic value is reached already for quite small wave-vectors is thus 
rather surprising. We note, however, that a closer inspection of _F s Na (g,t) for q > 1.6 A -1 
shows that the a-relaxation regime of these correlators fall into the late /3-relaxation regime 
of all the other (slow) correlators in which the von Schweidler law holds. Thus, it seems 
that the relaxation processes in the /3-regime of the slow correlators that correspond to 
the universal von Schweidler decay are impressed onto the dynamical behavior of the single 
particle motion of the sodium atoms which leads to stretched exponential decay of _F s Na (g, t) 
with (3 = b. For smaller wave-vectors this is not the case since there the a-relaxation of 
F s Na (g, t) overlaps essentially only with the plateau of the slow correlators. 

The result in Fig. ^ shows that the time dependence of the correlation functions for the 
time regime in which they start to fall below the plateau is compatible with the functional 
form given in Eq. @. The factorization property stated in Eq. @ is, however, more general, 
since expression (g) is just the leading asymptotic prediction for G(t). A different way to 
check to what extent the factorization property holds, without making use of the explicit 
form of G(t), is to calculate the following quantity: 

_ m - 

~ W) - W) ' { ] 

Here t' and t" are two arbitrary times in the /3-regime. It follows immediately that i/Eq. 
holds, is independent of </>, since it is just the system universal function G(t). In Fig. [7| 
we show the time dependence of at T = 2100 K. The correlators </> used are F s (q,t) for 
Si and O at q = 0.94 A -1 , q = 1.7 A -1 , q = 2.0 A -1 , and q = 3.0 A" 1 , as well as the 
coherent functions F(q,t) for Si-Si, Na-Na, and 0-0 at the same values of q. The times t" 
and t' from Eq. fllTf ) are 2.9 ps and 10.3 ps, respectively. From the figure we recognize that 
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in the /3-regime the R^t) for all these correlators collapse nicely onto a master function, 
which is the function G(t). Hence we conclude that the factorization property predicted by 
MCT holds for the present system. Finally we mention that we find for higher temperatures 
the same type of collapse, but that the time window in which the master curve is observed 
shrinks rapidly, in qualitative agreement with the prediction of MCT. 

The last quantity we will discuss is C from Eq. Since C is just the height of the 
plateau of the correlator at intermediate times, this parameter is often also called "non- 
ergodicity parameter" (NEP) since it reflects how much memory the system has of its state 
at t — 0. In the following we will focus on the wave- vector dependence of f al3 (q) and 
the NEP for the coherent and incoherent scattering functions. (Note that f a/3 {q) and ff{q) 
are often also called the Debye- Waller factor and Lamb-Mossbauer-factor, respectively.) 
We have determined f al3 (q) and ff(q) by using Eq. (|8|) to fit the correlators. The in- 
dependence of f a P{q) and ff{q) for the case of silicon and oxygen are shown in Fig. || We 
see that in both cases ff{q) (open symbols) decays quickly with increasing q and that this 
dependence can be described very well by a Gaussian (bold solid lines). Such a behavior 
has already been found for the case of pure silica and is in qualitative agreement with 
the prediction of MCT . The width q s of these Gaussians are 4.1 A" 1 and 3.0 A" 1 for Si 
and O, respectively. This means that in the time scale of the /3-relaxation the particles are 
trapped in a cage with radius r s = 1/q, = 0.24 A (Si) and 0.33 A (O). For the case of pure 
Si0 2 the corresponding values are 0.23 A and 0.29 A gf§. Hence 

we see that the presence of 

sodium slightly increases the size of the cage and that this increase is more pronounced for 
the case of oxygen than for silicon. This result is reasonable since some of the oxygen atoms 
are in the immediate vicinity of the sodium atoms but are bound to them less strongly than 
they are to the silicon atoms (and of course there is no Si-Na bond). In addition also the 
presence of dangling bonds (i.e. non-bridging oxygens) will lead to an increase of the size 
of the cage for oxygen. 

The NEP for the coherent functions of Si and O oscillate around the ones for the in- 
coherent functions. This oscillation is in phase with the corresponding structure factor, a 
behavior which is in qualitative agreement with the theoretical expectation || . In particular 
we see that the amplitude of this oscillation is smaller for the case of silicon than for oxygen, 



in agreement with the findings of pure Si02 |3(| or binary mixtures of particles [75 1. This 
finding can be rationalized by the fact that in a binary system with strong asymmetry in 
the concentration the coherent correlation functions for the minority species are in general 
very similar to the incoherent functions. 

Also included in the figure is f a (q) as determined from the Na correlator. (Note that 
F s Na (g, t) does not show a well defined plateau for q > 1.5 A -1 (see e.g. Fig. 2 in Ref. 
Hence a fit with the functional form given by Eq. (|8]) is rather difficult if q is large and thus 
the height of the plateau cannot be determined with high accuracy. However, for q < 1.2 A -1 
one does indeed find a well developed plateau and thus it is possible to determine with 
good accuracy.) From Fig. |8] we see that ff(q) for Na decays significantly faster than the 
ones for Si and O. This can be interpreted by saying that the cage for the sodium atoms is 
wider than the one for Si and O. However, care must be taken in drawing this conclusion 
since, as just mentioned, at the temperatures investigated neither the intermediate scattering 
functions for intermediate and large q nor the mean-squared displacement of the Na atoms 
show a well defined plateau at intermediate times. Hence one cannot really say that on 
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this time scale the particles are caged and hence also the interpretation of f£(q) as the 
Fourier-transform of the shape of the cage is not quite appropriate. What is remarkable 
with fg(q) for Na is the fact that it is not possible to fit it well with a Gaussian. If one fits 
the data for q > 2.2 A -1 with such a functional form it is possible to obtain a very good fit 
(bold solid line). (The width of this Gaussian is 2.52 A -1 , which corresponds to a "cage" 
of size 0.39 A.) However, this fit gives a very poor representation of the data for smaller 
wave-vectors. This result, which is in contrast to the findings for Si and O, shows that the 
dynamics of Na is indeed rather unusual. Roughly speaking one thus can say that for small 
length scales, i.e. large q, the cage is relatively soft and therefore the ff a (q) decays quickly. 
However, on the length scales of the typical distances between the channels (small q) the 
"cage" is relatively rigid, since on this length scale, as we have mentioned before, the single 
particle dynamics of Na is not strongly coupled to relaxation processes in the matrix. But 
of course this is just a hand-waving explanation of the finding and it would be nice to find 
support for it by means of a more thorough theoretical calculation. 

Also included in the figure is the NEP for the Na-Na correlation. Qualitatively this 
f a/3 (q) looks similar to the one for the Si-Si or 0-0 correlation and in particular it shows 
a pronounced peak at around 2.0 A" 1 , which corresponds to the nearest neighbor distance 
between two Na atoms (= 3.3 A) |]66|| . However, we notice an important difference in that 
this NEP does not oscillate around the NEP for the incoherent function but instead stays 
systematically below it. Also for this behavior we are not aware of any theoretical prediction 
or experimental result. 



IV. SUMMARY 

In this paper we have studied by means of molecular dynamics computer simulations the 
relaxation dynamics of a melt of Na20-2Si02, one of the prototypes of an ion-conducting 
glass former. In particular we investigated the temperature and wave-vector dependence 
of the coherent and incoherent scattering functions. Due to the very different time scale 
of the dynamics of Na from the one of the species forming the matrix, the g-dependence 
of these correlators for the sodium atoms shows features that are neither found in simple 
liquids nor in network-forming liquids like pure silica, i.e. systems in which the dynamics 
of the individual species takes place on a comparable time scale. E.g. we find that the time 
and temperature dependence of the incoherent function for Na is very different from the one 
for the coherent function. This is related to the fact that in this system the sodium atoms 
have the tendency to populate a relatively small subregion of space, so-called channels, 
and that the dynamics of the atoms in these channels is relatively quick and occurs by 
(activated) single particle hops. In contrast to this the overall structure of the channel, 
and hence the coherent function for the Na atoms, relaxes only on the time scale of the 
a-relaxation time of the Si-0 matrix. The temperature dependence of the relaxation times 
for the matrix and _p NaNa (g ; t) show a strong deviation from an Arrhenius law in agreement 
with experimental findings for this system [|79| , p0[ . In the past such deviations have been 



found for simple liquids [|75l , pIH 83[| and their existence has been rationalized by means of 



MCT. However, finding them in systems like silica [^,^] or in the present sodium silicate 
system is rather surprising. In this paper we have shown that in the temperature regime 
where these deviations are seen many of the features of the relaxation dynamics can again be 
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rationalized by means of MCT, thus showing that with respect to this there is no difference 
with the results found for the so-called fragile glass-formers. Hence we conclude that the 
main difference between strong and fragile glass-formers is that the presence of the hopping 
processes leads to a shrinking of the dynamical range in which the T-dependence of the 
relaxation times follows the power-law predicted by MCT (as compared to the range found 
in fragile systems). Nevertheless, despite this reduced range, the time correlation functions 
still show a behavior that can be rationalized remarkably well by the theory. 

Of course one has to wonder to what extent the results presented in this paper can be 
found also in real Na20-2Si02, or similar systems. Although it must be expected that the 
potential used is not sufficiently accurate to reproduce all the properties of the real material 
on a quantitative level, the surprisingly good agreement of the results of the present model 



with the neutron-scattering results of Meyer et al. [[M] shows that the potential is quite 
realistic. Therefore it can be hoped that the results presented here will be found also in 
appropriate coherent and incoherent neutron-scattering experiments and we hope that the 
present work helps to motivate such experiments. 
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FIG. 1. Time dependence of the incoherent intermediate'scattering function for Si, O, and Na 
(top to bottom) for all temperatures investigated and a wave-vector q = qi = 1.7 A -1 . The dashed 
lines are fits to the curve for T = 2100 K with a KWW-law. 
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FIG. 2. Time dependence of the coherent intermediate scattering function for the Si-Si, (a), 
and Na-Na, (b), correlation for all temperatures investigated. 
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FIG. 3. The incoherent function F s (q,t) versus t/r(T) for Si, Na, and O (top to bottom). The 
a-relaxation time r(T) has been defined via F s (q, r) = 0.1. The inset in (a) and (b) show the same 
functions at large rescaled times. 
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FIG. 4. Wave-vector dependence of rq 2 for various coherent and incoherent correlators (see 
labels of the curves) at T = 2100 K. Right scale: g-dependence of the partial structure factor for 
the sodium atoms. 



18 




1800 2400 3000 3600 4200 

T[KJ 

FIG. 5. Test for the presence of a power-law for the relaxation times r{q) for Si and O and 
two wave-vectors. Also included is the data for the diffusion constants Dgi and -Do, which for the 
sake of a clearer presentation have been multiplied by 35. The bold straight lines are linear fits 
corresponding to the theoretical expectation. 




FIG. 6. Time dependence of the incoherent intermediate scattering function for oxygen at 
T = 2100 K. The solid curves are F®(q,t) for different wave-vectors (see labels). The dotted 
curves are fits with a von Schweidler law, first term in Eq. (|9|), with an exponent b = 0.47. The 
dashed curves are fits with the von Schweidler law including the leading order corrections (see 
Eq. (I)). 
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FIG. 7. Time dependence of the function R<j,(t), denned in Eq. QlOj) , at T = 2100 K. The 
correlators used to make this plot are Fg*(q, t) and F®(q,t) for q = 0.94,1.7,2.0,3.0 A -1 and 
F aa (q,t) with a € {Si, 0,Na} and at the same wave-vectors. 
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FIG. 8. Wave-vector dependence of the the non-ergodicity parameters f aa (q) and ff{q) 
(curves with open and closed symbols, respectively) for a € {Si, 0,Na}. The bold lines are Gaus- 
sian fits to the NEP for the incoherent functions. 
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